plot_color

r = GET('http://data.beta.nyc//dataset/3bf5fb73-edb5-4b05-bb29-7c95f4a727fc/resource/6df127b1-6d04-4bb7-b983-07402a2c3f90/download/f4129d9aa6dd4281bc98d0f701629b76nyczipcodetabulationareas.geojson')
nyc_zipcode = readOGR(content(r,'text'), 'OGRGeoJSON', verbose = F)

## No encoding supplied: defaulting to UTF-8.

#summary(nyc_zipcode)

#nyc_zipcode_df <-= tidy(nyc_zipcode)

leaflet(nyc_zipcode) %>%
  addTiles() %>% 
<<<<<<< HEAD
  addPolygons(popup = ~neighborhood) %>%
  addProviderTiles("CartoDB.Positron")

## add points
=======
  addPolygons(popup = ~postalCode)

geo_jion

data_hiv = read.csv("./data/combine_hiv_new.csv") %>%
  filter(gender == "All", age == "All", race == "All")

diagnoses_by_zipcode = data_hiv %>% 
  group_by(zipcode, uhf) %>% 
  summarise(sum_hiv_diagnoses = sum(hiv_diagnoses))

map_data = geo_join(nyc_zipcode, diagnoses_by_zipcode, "postalCode", "zipcode")

point

#nyc_zipcode_df = as_tibble(nyc_zipcode@data) 
#points_spdf = full_join(data_hiv, nyc_zipcode_df, by = "zipcode") %>% 
>>>>>>> fd8234d54b79076b674a1f6a542c09da1cda0497

points_spdf = read.csv("./data/combine_hiv_pin.csv")

coordinates(points_spdf) = ~lng + lat
proj4string(points_spdf) = proj4string(nyc_zipcode)
matches = over(points_spdf, nyc_zipcode)

diagnoses_by_postalcode = diagnoses_by_zipcode %>% 
  rename(postalCode = zipcode)

points = full_join(matches, diagnoses_by_postalcode, by = "postalCode")

## Warning: Column `postalCode` joining factors with different levels,
## coercing to character vector

color

pal = colorNumeric(palette = "plasma", domain = range(map_data@data$sum_hiv_diagnoses, na.rm = T))
# "BuPu" "viridis" "Greens""inferno"

leaflet(map_data) %>%
  addTiles() %>%  
  addPolygons(color = "black", weight = 0.5, opacity = 1, 
              fillColor = ~pal(sum_hiv_diagnoses), 
              popup = ~stringr::str_c(uhf, "  sum:", factor(sum_hiv_diagnoses))) %>% 
  addMarkers(~longitude, ~latitude, 
             popup = ~stringr::str_c(uhf, "  sum:", factor(sum_hiv_diagnoses)), data = points) %>%  
  addProviderTiles("CartoDB.Positron") %>%
<<<<<<< HEAD
  setView(-73.98, 40.75, zoom = 13)

======= setView(-73.98, 40.75, zoom = 11)

## Warning in validateCoords(lng, lat, funcName): Data contains 124 rows with
## either missing or invalid lat/lon values and will be ignored

>>>>>>> fd8234d54b79076b674a1f6a542c09da1cda0497